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Abstract. Motivated by the rapidly growing possibilities for experiments with 
ultracold atoms in optical lattices we investigate the thermodynamic properties of 
correlated lattice fermions in the presence of an external spin-dependent random 
potential. The corresponding model, a Hubbard model with spin-dependent local 
random potentials, is solved within dynamical mean-field theory. This allows us to 
present a comprehensive picture of the thermodynamic properties of this system. In 
particular, we show that for a fixed total number of fermions spin-dependent disorder 
induces a magnetic polarization. The magnetic response of the polarized system differs 
from that of a system with conventional disorder. 
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1. Introduction 

Interacting quantum many-particle systems in the presence of disorder are of great 
interest not only in condensed matter physics [H [21 El IH O |6] but recently also in 
the field of cold atoms in optical lattices [7l[8l|9l[ini[IIl[l2l[l3l[14]. In particular, 
ultracold gases have quickly developed into a fascinating new laboratory for the study 
of quantum many-body physics. A major advantage of cold atoms in optical lattices is 
the high degree of controlability of experimental parameters such as the interaction and 
the disorder strength, which allows for a detailed comparison between experiment and 
theory. Indeed, advanced techniques developed for experiments with ultracold gases in 
optical lattices facilitate the realization and explicit study of model situations which are 
not accessible in solid state physics. For example, it was recently demonstrated that 
it is possible to prepare optical lattices with a spin- dependent hopping of the quantum 
particles [151 IIHI [H]. Corresponding theoretical models had already been proposed 
earlier [HI UHl EQl EH ESI ES] . 

In this paper we discuss the thermodynamic properties of interacting spin-1/2 
lattice fermions in the presence of a spin-dependent random potential. Such a model 
was recently investigated by Scalettar and collaborators [21] . It may be realized 
experimentally in the following way: In an off-resonant electric field atoms experience 
Stark shift of their ground state energies, which is proportional to the light intensity. 
The optical lattice is formed by superimposing laser waves propagating in different 
spatial directions. This effectively leads to a periodic dependence of the electric field 
intensity and thereby to a periodic potential. Disorder can be introduced by focusing 
a light beam on the optical lattice, which is scattered from a diffusive plate. Such 
an effective random potential (speckle-type disorder) is characterized by a very short 
correlation length and a pronounced statistical independence [HI [121 [25]. Using laser 
beams with different polarization it is possible to generate a random potential which 
acts differently on particles with different spin-orientation in the ground state |24j . 
The spin-dependent speckle disorder affects the local one-particle energies, hopping 
amplitudes, and interparticle interaction potentials. Since ultracold atoms in optical 
lattices can be experimentally prepared in quite different ways, two model situations 
are of particular interest: (a) the total number of fermions is conserved, and (b) the 
number of fermions with different spin projection is individually conserved; the latter 
case is refered to as a "spin-imbalanced" system. Nanguneri et al. [21] analyzed case (b) 
for an attractive interaction between the fermions, which leads to a pairing instability. 
Employing Bogoliubov-de Gennes mean-field theory they found that spin-dependent 
disorder leads to a suppression of the pairing state. 

In the first part of this paper we investigate the case when the total number of 
fermions is conserved (i.e., when the two spin populations have a common chemical 
potential) for a repulsive local interaction. In particular, we will focus on fillings off 
half-filling, since in this case spin- dependent disorder induces a finite magnetization 
through the broadening of the spin sub-band on which the disorder acts. At the same 
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time characteristic correlation features which occur only at half-filling, such as the Mott 
metal- insulator transition [H |26l [271 [2H1 l29l |30] , will be absent |3T]. The spin- dependent 
randomness is introduced through local, random, and spin-dependent potentials. Such 
a model is closely related to electronic models investigated in solid state physics [31], 
in particular correlated lattice electrons in random potentials; see, for example, refs. 
[35| |36| [371 IMl |39l HQ]. However, in sohd state systems it is not possible (or, at least, 
not yet possible) to tune and manipulate spin- dependent potentials experimentally. 
By contrast, as discussed above, optical lattices with ultracold atoms do offer such 
possibilities. In particular, the model where the total number of fermions is conserved 
is applicable to systems with atoms characterized by a fast dipolar relaxation, such as 
Cr or rare elements with large magnetic dipole moments [41j. Similarly, by employing 
double-shot phase contrast imaging for atoms in an external parabolic potential trap 
it is possible to determine spin- dependent densities, the spin susceptibility, and the 
compressibility of particles [12]. Hence systematic studies of the disorder induced 
spin-asymmetric changes of the thermodynamic properties of such system are possible. 
In the second part of the paper we report results for spin-imbalanced systems, i.e., 
when each spin population has an individual chemical potential. Spin-imbalanced 
fermionic systems have been investigated experimentally |^ [JH [15] and theoretically 
[SI HZl HH Hni [50] either in optical lattices or in the case of a trap without disorder. 
Spin-dependent disorder adds a new tool to the examination of correlation effects in 
spin-imbalanced fermionic systems. 

2. Hubbard model and spin-dependent disorder 

In the following we will study interacting fermions in a random environment, modeled 
by the Anderson-Hubbard Hamiltonian with local, spin-dependent disorder 

ij,(7 ia icr i 

Here Co- (c|^) is the anihilation (creation) fermionic operator, hia^ = cl^Ci^j is the particle 
number operator, tij is the hopping matrix element, U is the local interaction, and fia is 
the chemical potential for particles with spin a. In the model where the total number of 
fermions is conserved the chemical potentials of the two species obey fi^ = = /i. In 
the model where the numbers of fermions with different spin are individually conserved 
the chemical potentials fi„ can differ. 

The disorder is represented by random local potentials which are distributed 
according to a box probability density P{eia) = l/A^ for |ejo-| < ^0-/2, and zero 
otherwise ("continuous disorder"). Here Ao- is the maximal energy difference between 
the local energies for a given spin direction and thus provides a measure of the disorder 
strength. Averages over the disorder are calculated by (■ ■ ■)dis = 110-/ dea^P{e^){- ■ ■), 
whereby the study of Anderson localization within the one-particle Green function 
formalism is excluded [511 [521 [53] . 
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It is interesting to note [m [Mj that by introducing the z-component of the 
magnetization operator rhi = fii^ — hi^ and the local particle number operator hi = 
^it + Anderson-Hubbard Hamiltonian for the model where the total number of 

fermions is conserved takes the form 

H = J2 tijcLcja - J2 Ih^i + K^i + fli^fli^, (2) 

ij,a i i i 

where jj,* = ji — {ei^ + eii)/2 and h* = (cj-f- — eii)/2. The Hamiltonian ([2]) has a natural 
interpretation: it describes fermionic atoms moving in a random chemical potential /x* 
and a random Zeeman magnetic field h* which are correlated, since 

(/x*/.*) = ^(A5-AX-. (3) 

When Ai = A-]- we arrive at the Anderson- Hubbard model with spin-independent 
disorder [55] . 



3. Dynamical mean-field theory 

We solve the Hamiltonian within dynamical mean-field theory (DMFT) |56l 
Hence all local dynamical correlations due to the local interaction are fully taken into 
account. However, non-local correlations in position space are absent. Long-range order 
in space can be included within the DMFT, but this will not be the case in the present 
investigation. 

In the DMFT scheme the local Green function Gan is determined by the bare density 
of states (DOS) N^{e) and the local self-energy Eo-„ as Gan = J deN^{e) / {iun + fia — 
^an — e)- Here the subscript n refers to the Matsubara frequency iun = ii^n + l)7r//3 
for the temperature T, with /3 = 1/T. Within DMFT the local Green function Go^n is 
determined self-consistently by 

ri — —llr r* W — J ^ i^^'^ ^cr\ ^an^an^ \ /.n 

together with the k-integrated Dyson equation Q^^ = G^^ + S^n ^281 157]. The single-site 
action Ai for a site with the ionic energy G [— 4; t] has the form 



★ 1 * 1 f'^ 

n,u a 

-^E [' dTc:{T)cAr)c%{r)c_^{r), (5) 
2 ^ Jo 

where we used a mixed time/frequency representation for Grassmann variables c„, c* . In 
Eq. mthe average over the quantum states with Boltzmann distribution (inner bracket) 
and over the disorder (outer bracket) is taken. In the non-interacting case [U = 0) the 
DMFT equations reduce to those of the coherent potential approximation [581 Ell |60] . 



In the following we employ a semicircular DOS A^°(e) = {8/'7iW^)yW'^/4 — e'^, 
which corresponds to nearest-neighbor hopping on a Bethe lattice with infinite 
coordination number (Z = oo) and choose the half band width W = 2 as the 
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energy scale. The one-particle Green function in Eq. is determined by solving the 
DMFT equations iteratively [611 162] using Hirsch-Fye Quantum Monte-Carlo (QMC) 
simulations [63] . 

We now discuss and compare results for different types of disorder: i) spin- 
independent disorder, where = = A, and ii) spin- dependent disorder with A^ = A 
and A-i- = 0. In both cases A is a measure of the disorder strength. They represent 
extreme situations and thereby illustrate the most important differences between spin- 
independent and spin- dependent disorder. Intermediate cases can be studied in the 
same way. 

In the model where the total number of fermions is conserved the chemical potential 
fi is fixed to keep the total number of fermions per site constant: n = + n^, where 
^cr = fiia))dis/NL and Nl is the number of lattice sites. The numerical investigations 
were performed for the densities n = 0.3, 0.5 and 0.7 pi]. Qualitatively, all three cases 
lead to similar results. For the purpose of illustration we present results only for n = 0.5. 

In the model where the numbers of particles with different spin a are separately 
conserved the chemical potentials /i^ are fixed individually to keep the number of 
fermions with spin a per site constant. The numerical investigations were performed for 
(n^,nj = (0.2,0.3) , (0.3,0.2), (0.2,0.5), and (0.5,0.2). For the purpose of illustration 
we here present results for (0.2,0.5), and (0.5,0.2) since qualitatively, all cases lead to 
similar thermodynamic properties. 

4. Constant total number of fermions 

4.I. Finite magnetization caused by spin- dependent disorder 

The most striking effect caused by spin-dependent disorder is the appearance of a finite 
magnetization m(A, T, U) = '^i))dis/^L for A > at arbitrary temperatures T and 
interaction values U, as shown in the upper panel of Fig. [T] for the lowest temperature 
studied here, T = 0.06. The magnetization m is seen to grow monotonously as the spin- 
dependent disorder strength A is increased. The increase of m is slightly stronger for 
larger values of the interaction U. It should be noted that the origin of the magnetization 
is not due to electronic correlations, but is a pure one-particle effect, which for a 
symmetric DOS, occurs only away from half filling (n 7^ 1). Namely, spin-dependent 
disorder leads to a symmetric broadening of the spin subband on which it acts (here the 
down-spins), whereby the number of spin-down particles increases while the number of 
spin-up particles decreases since the total number of particles remains constant (that is, 
the two spin populations have a common chemical potential as shown in the lower part 
of Fig. [1]). This then leads to a finite magnetization. Here we assumed that the system 
is able to relax to the new equilibrium state by means of spin-flip processes which are 
provided, for example, by the spin-orbit interaction or inelastic dipolar collisions |41j . 
but which are not explicitly included in the model Hamiltonian ([T]). The weak influence 
of U on the magnetization may be attributed to the flow of spectral weight induced by 



Correlat 




Figure 1. Conserved total number of fermions: Upper panel: Increase of the 
magnetization m as a function of spin-dependent disorder A for interactions U = 
0.1,1.0,4.0 at T = 0.06. Similar results are obtained from a non-interacting model 
where the band widths of particles with opposite spin differ by an energy S (see text). 
Lower panel: Schematic drawing of the influence of spin-dependent disorder. The 
disorder leads to a broadening of the respective spin subband. Due to the conservation 
of the total number of fermions a finite magnetization is built up. 

the correlations. 

For a symmetric DOS and a half-filled band the number of states below the Fermi 
energy Ep = remains constant when A is increased; this holds for any interaction 
U. The magnetization is then zero. This effect was clearly observed in our numerical 
simulations, but is not presented here. 

The increase of the magnetization m with increasing spin- dependent disorder A 
may be understood already within a simple model without interactions {U = 0) where 
the particles with different spin orientation have different band widths [65l |66l |67] . 
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Namely, if we assume that the DOS of the system is given by 

f J_ for \f\ < ^ 
'^^'^ \ for \e\>^, 
where = W and = W + 6j the magnetization at T = is given by 

m(n,5) = (l-n)^^^. (7) 

The increase of the magnetization as a function of the difference of the band widths, 
(5, at T = is shown by the full black curve in the upper panel of Fig. [H A very 
similar result, indicated by a black dashed curve in the same figure, is obtained for 
a semi-elliptic DOS with spin-dependent band widths N^{e) = {S/TrW^)^W^/4: — e^. 
Both curves closely resemble those calculated for the Hamiltonian ([1]), except for the 
turning point observed in m{A,T,U) for U = 0.1,1.0. Apparently this feature is due 
to scattering effects on the impurities which are included in the self-energy of the local 
Green function. 



4-2. Thermodynamic properties 

As the temperature increases the Fermi-Dirac distribution is smeared out and the 
arguments for the appearance of a finite magnetization due to a spin-dependent 
band broadening at a constant chemical potential, although still valid, become less 
stringent. Indeed, as shown in Fig. [2]the magnetization slowly decreases with increasing 
temperature for all values of the interaction U and disorder strengths A > 0. 



Correlated fermions in lattices with spin- dependent disorder 



8 



The magnetization induced by the disorder influences the thermodynamic 
properties of the system. In Fig. |3] we compare the temperature dependences of the 
average double occupation d = dis/ Nl for spin-independent disorder (left 

panel) and spin-dependent disorder (right panel), respectively. For spin-independent 
disorder and weak interactions, e.g. U = 0.1, an increase of the disorder A leads 
to a corresponding increase of d, i.e., on average more energy levels below the Fermi 
energy are occupied by two fermions with opposite spins. For increasing U a minimum 
appears in d{T) which corresponds to a maximum in the local moment per site 
S = ^J {J2i{'>T^'i))dis/NL = y/n — 2d. This is seen clearly for f/ = 4 in Fig. [31 Increasing 
the disorder A reduces the value of the local moment even in the presence of strong 
interactions (f/ = 4), when d is very small. Only at very low temperatures and U = 4 
we can see an opposite trend. 

In the case of spin-dependent disorder (right panel of Fig. [3]) the double occupation 
is always found to be suppressed by the disorder A. This is an effect of the finite 
magnetization of the system. As a consequence the local moment increases. We also 
see that the local minimum of d is shifted to lower temperatures when A is increased. 
Whether it disappears at T = cannot be decided here since we are not able to run 
QMC simulations at lower temperatures. 
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Figure 3. Constant total number of fermions: Double occupation d as a function 
of the temperature T for different values (A = 0, 1, 3, 5) of spin-independent disorder 
(left) and spin-dependent disorder (right); upper panel: U = 0.1, middle panel: U = 1, 
lower panel: U — A. 
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Figure 4. Constant total number of fermions: Density of states at the chemical 
potential as a function of the temperature T for different values (A = 0, 1, 3, 5) of spin- 
independent disorder (left) and spin-dependent disorder (right); upper panel: U = 0.1, 
middle panel: C/ = 1, lower panel: U = A. 



Next we compare the value of the DOS at the chemical potential N{fi) for different 
disorder strengths (Fig.H]). As expected, an increase of the disorder reduces the DOS. We 
also see that this decrease is stronger for spin-independent disorder because the random 
potential acts equally on both spin subbands. When the temperature is lowered, the 
DOS increases similarly in both cases due to the enhanced quantum coherence of the 
system. 

Static susceptibilities provide useful information about possible phase instabilities 
of the system and their response to external fields. Therefore, we now compute 



the ferromagnetic (FM) susceptibility Xfm ~ '^^ 



dh 



, where m is the magnetization 



dnis 



where m 



St 



density in the presence of an external magnetic field h 
the antiferromagnetic (AFM) susceptibility Xafm — 

riA^ — is the staggered magnetization density on a bipartite lattice with 
nonequivalent sites A and B in the presence of a staggered magnetic field hst, 
where nA{B) = {{EieA{B) '"'ia)) dis, and 

the density susceptibility (compressibility) Xc = where n is the particle number 
density at the chemical potential fi. 
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Figure 5. Constant total number of fermions: Inverse ferromagnetic susceptibility as 
a function of the temperature T for different values (A = 0, 1, 3, 5) of spin-independent 
disorder (left) and spin-dependent disorder (right); upper panel: U = 0.1, middle 
panel: ?7 = 1, lower panel: C/ = 4. 



The magnetic response in the FM and AFM channels differ significantly for spin- 
independent and spin-dependent disorder, respectively: both susceptibilities are reduced 
by spin-independent disorder, at least for weak interactions [U = 0.1, 1), as seen in the 
upper and middle left panels of Figs. |5] and |6l For strong interactions in the presence of 
spin-independent disorder {U = 4, lower left panels of Figs. [5]and|6]) the FM and AFM 
susceptibilities behave differently, in particular at the lower temperatures. Namely, 
they exhibit Curie-like behavior and increase with increasing disorder. This is due 
to the formation of local moments and corroborates the temperature behavior of the 
double occupation discussed earlier. Indeed, at low enough temperatures and for strong 
interactions the double occupation d is reduced by spin-independent disorder. 

A similar Curie-like behavior of the FM susceptibility was previously found in the 
two-dimensional Hubbard model with random box potential [68]. However, in this 
case neither the AFM susceptibility nor the averaged local moment show any unusual 
behavior. Therefore we conclude that in contrast to high dimensions, where correlated 
electrons in random potentials form local moments, the behavior in two dimensions is 
different and is yet unexplained. 

In the case of spin-dependent disorder an increase of the disorder strength A reduces 
both magnetic susceptibihties (see the right panels of Figs. [5]and|6]). Namely, due to 
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Figure 6. Constant total number of fermions: Inverse antiferromagnetic susceptibility 
as a function of the temperature T for different values (A = 0, 1, 3, 5) of spin- 
independent disorder (left) and spin-dependent disorder (right); upper panel: U = 0.1, 
middle panel: U — \, lower panel: C/ = 4. 



the magnetic polarization induced by the spin- dependent disorder the response of the 
system to a magnetic field is weaker. 

In Fig.[7]we compare the temperature dependence of the compressibility of a system 
in the presence of spin-independent and spin- dependent disorder, respectively. In both 
cases disorder reduces the compressibility and, as in the case of the DOS, the effect of 
spin-independent disorder is stronger. We also observe an enhancement of the statistical 
fiuctuations within the Monte Carlo method at stronger interactions. 

Finally we note that at any given temperature T the reduction of the compressibility 
Xc by disorder is smaller in the strongly correlated case than in the weakly correlated 
case. This holds true for both types of disorder, and is due to the fact that a strong 
repulsive interaction leads to a rigidity of the system which makes it less sensitive to 
the infiuence of disorder. 

The non-interacting model discussed in Section [3] where particles with different spin 
have different bandwidths [65] yields qualitatively similar results for the ferromagnetic 
susceptibility and the compressibility. Although such a toy model cannot explain the 
details of the magnetization it appears to be quite useful for explaining qualitative 
features of systems with spin-dependent disorder. 
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Figure 7. Constant total number of fermions: Compressibility Xc as a function of the 
temperature T for different values (A = 0, 1,3, 5) of spin- independent disorder (left) 
and spin-dependent disorder (right); upper panel: U = 0.1, middle panel: U = 1, lower 
panel: U ~ A. 



5. Spin-imbalanced fermions 

In the case of spin-imbalanced fermions the number of fermions with spin a is conserved 
individually. The magnetization m = — is then constant, i.e., does not 
depend on the thermodynamic variables and parameters of the model. In Fig. |8] we 
compare the temperature dependence of the average double occupation d of the spin- 
imbalanced fermions for spin-independent disorder (left column) and spin-dependent 
disorder (middle and right columns), respectively. For spin-independent disorder in the 
presence of a weak interaction {U = 0.1) an increase of the disorder strength A leads to 
a corresponding increase of d, i.e., on average more energy levels below the Fermi energy 
are occupied by two fermions with opposite spins. For increasing U a minimum appears 
in d{T) which corresponds to a maximum in the local moment. This is seen clearly for 
[/ = 4 in the left column of Fig. [HI Therefore a maximum of the of local moments is 
found in both models considered here. In the case of spin-dependent disorder (middle 
and right columns of Fig. [8]) the double occupation is, in general, suppressed by the 
disorder A. This is an effect of the finite magnetization of the system, which is here 
imposed by fixing the individual spin densities. As a consequence the local moment 
increases. We also see that the local minimum of d{T) is shifted to lower temperatures 
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Figure 8. Spin-imbalanced fermions: Double occupation d as a function of the 
temperature T for different values of A — 0,1,3,5 for spin- independent disorder 
(left column) and spin-dependent disorder (middle and right columns); upper panels: 
U = 0.1, middle panels: U — 1, lower panels: ?7 = 4. In the left and the middle 
columns the density of particles is — 0.2 and — 0.5 and in the right column the 
density of particles is = 0.5 and = 0.2. 



when A is increased. Whether it disappears at T = cannot be decided here since, as 
mentioned earher, we are not able to run QMC simulations at lower temperatures. It is 
interesting to note that the behavior of the double occupation is qualitatively similar in 
both models. 

In Fig. [9] we compare the values of the spin resolved DOS at the chemical potential 
for different types of disorder and different disorder strengths. Spin-independent disorder 
always reduces the DOS, while spin-dependent disorder only reduces the DOS of the 
corresponding spin subsystem. For weak interactions the DOS of the opposite spin 
subsystem is almost unchanged. Only at stronger interactions {U = 1 and f/ = 4) 
does the disorder in one spin subsystem also influence the DOS of the opposite spin 
subsystem (but only weakly). This correlation effect is expected, since the two spin 
subsystems are coupled by the on-site interaction U . 

In the case of spin-imbalanced fermions it is useful to discuss a compressibility 
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Figure 9. Spin-inibalanced fermions: Spin resolved density of states for spin down 
(upper figure) and spin up (lower figure) at the chemical potentials as a function of 
the temperature T for different values (A = 0, 1, 3, 5) of the spin-independent disorder 
disorder (left) and the spin-dependent disorder (right); upper panels: U = 0.1, middle 
panels: U = 1, lower panels: U — A. In the left and the middle columns the density of 
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matrix [SS] 



Xc 




T 



(8) 



where n^. = X^i^io-- Its off-diagonal elements provide a measure of density-density 
correlations between different spin subsystems. We note that the matrix is symmetric 
Xcti = Xcit- AtU = one has Xcu = 0- 

In Fig. [To] the diagonal compressibilities in different spin channels are shown. In the 
case of spin-independent disorder (left columns in Fig. [TOj) the diagonal compressibilities 
are found to be reduced in both spin subbands for increasing disorder strength A. Due 
to the imbalance of the spin population one has Xctt XcU- 

In the case of spin- dependent disorder the spin-down compressibility Xc u reduced 
when the disorder A is increased; this holds for any interaction strength U, cf. upper 
figure, middle and right columns in Fig. [TOl By contrast, the spin- up compressibility 
Xctt remains almost unchanged at weak interactions. But for larger U it decreases 
for increasing disorder A, cf. lower figure, middle and right columns in Fig. [TOl The 
disorder-induced changes in Xctt interaction effect, since the disorder A acts only 

on the spin-down subsystem. 

For spin-imbalanced fermions we observe a distinctive qualitative difference between 
spin-independent and spin-dependent disorder. Namely, as for the model where only the 
total number of fermions is kept constant, we find that at any given temperature T the 
reduction of the compressibility Xc by spin-independent disorder is less in the strongly 
correlated case than in the weakly correlated case (quantitatively this effect is even 
much stronger here, cf. Fig. [7]). But the opposite happens in the case of spin-dependent 
disorder: then the reduction of the compressibility by disorder increases with increasing 
interaction strength U. 

In Fig. [H] we show the off-diagonal elements of the compressibility. As expected 
their absolute values increase with increasing interaction U. We also observe that for 
both types of disorder an increase of the disorder strength leads to a reduction of 
the inter-spin correlations. Although the off-diagonal elements of the compressibility 
are negative the determinant of the compressibility matrix is positive. Therefore, 
the system is stable against phase separation. The effect of disorder on the off- 
diagonal compressibilities at a given temperature T is similar to that of the diagonal 
compressibilities. 

6. Conclusions and Outlook 

In summary, we explored the thermodynamic properties of correlated fermions in the 
presence of spin-independent and spin- dependent disorder within dynamical mean-field 
theory. We discussed two models where either the total number of fermions, or the 
number of each spin component, is conserved. These two cases can be realized in 
experiments on cold atoms in optical lattices. 
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Figure 10. Spin-imbalanced fermions: Diagonal compressibilities Xc (upper figure) 
and Xc tt (lower figure) as a function of the temperature T for different values 
(A = 0, 1, 3, 5) of spin-independent disorder (left columns) and spin-dependent disorder 
(middle and right columns); upper panels: U — 0.1, middle panels: U — 1, lower 
panels: U — A. In the left and the middle columns the density of particles is = 0.2 
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Figure 11. Spin-imbalanced fermions: Off-diagonal compressibilities Xc it ^ 
function of the temperature T for different values (A = 0, 1, 3, 5) of spin- independent 
disorder (left columns) and spin-dependent disorder (middle and right columns) ; upper 
panels: U = 0.1, middle panels: U = 1, lower panels: U = 4. In the left and the middle 
columns the density of particles is = 0.2 and n-f- = 0.5 and in the right column the 
density of particles is = 0.5 and Uf = 0.2. Note that the relative error bars are 
of the same order in all cases but the off-diagonal susceptibility is very small at weak 
interactions, cf. compare scales on vertical axis. 

In the first model we found that, in contrast to spin-independent disorder, spin- 
dependent disorder induces a ferromagnetic polarization. However, instead of a Zeeman- 
Stoner-type shift of the sub-bands a spin-dependent broadening results. 

In the second model we showed that disorder which acts only on fermions with 
one particular spin direction nevertheless also affects the properties of particles with 
the opposite spin direction; this effect is more easily observed at strong interactions. 
Therefore, the investigation of the spin-resolved densities of states and spin-resolved 
compressibilities will provide useful quantitative information about correlations between 
the two spin subsystems. 
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